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Abstract. We demonstrate interesting and previously unforeseen properties of a pair of gap solitons in a 
resonant photonic crystal which are predicted and explained in a physically transparent form using both 
analytical and numerical methods. The most important result is the fact that an oscillating gap soliton 
created by the presence of a localized population inversion inside the crystal can be manipulated by means 
of a proper choice of bit rate, phase and amplitude modulation. Developing this idea, we are able to 
obtain qualitatively different regimes of a resonant photonic crystal operation. In particular, a noteworthy 
observation is that both the delay time and amplitude difference must exceed a certain level to ensure 
effective control over the soliton dynamics. 

PACS. 42.70.Qs Photonic bandgap materials - 42.65.Tg Optical solitons; nonlinear guided waves - 
42.50.Md Optical transient phenomena - 03.65.Ge Solutions of wave equations: bound states 



1 Introduction 

There has been an extensive amount of research focused 
on the study on nonlinear periodic systems and in partic- 
ular those described by the sine-Gordon equation. This is 
primarily due to a large and diverse range of physical phe- 
nomena which can be described by the soliton solutions of 
this equation. Resonance photonic crystals (RPhC) rep- 
resent one such system where the sine-Gordon equation 
can be used to describe the propagation of light. Periodic- 
ity in photonics results in bandgaps, or regions in time or 
space, where the propagation of energy is forbidden [TJ[U 
[5]. The addition of nonlinearity to the system allows the 
possibility of locally detuning these band gaps with the 
result that energy can propagate in the form of a solitary 
wave, or gap soliton (GS) [H|5j[6j[7] . Such gap solitons have 
been experimentally observed in fiber Bragg gratings [H] 
and more recently in AlGaAs waveguides [9 UTDllll| . These 
experiments relied on the non-resonant, Kerr nonlinear- 
ity to locally detune the stop-bands associated with the 
periodic structure. 

A recent resurgence of interest in photonics crystals 
with resonant nonlincaritics is due to the advances in the 
manufacturing of periodic Bragg semiconductor nanos- 
tructures such as InGaAs/GaAs multiple quantum wells 
and rare-earth doped AlGaAs/GaAs structures [T^llTSllbl] . 
and understanding their linear and nonlinear optical prop- 
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erties [T31HDl[r7] . This suggests that the building blocks of 
future all-optical processing systems, such as those based 
on the propagation of GSs might be experimentally feasi- 
ble at moderate optical power levels of about 10 MW/cm 2 . 
The inclusion of a gain or a local population inversion in 
the system opens up a range of additional, novel nonlinear 
effects PI]. 

Recently, we proposed an effective method for the con- 
trol of speed of gap solitons based on the defect associated 
with a weak linear excitation, or incoherent pump inside 
an RPhC [Tgil3fMn i23] . It turns out that a break in pe- 
riodicity of the RPhC is not necessary for slowing a gap 
soliton down and its trapping because the localized field 
creates a potential that acts as a phase-sensitive trap. The 
GS interaction at such defect was first studied to deter- 
mine the gain implication for slow light localization in the 
RPhC pjj| . More recent work was focused on the effects of 
localized gain on the GS collision and comprehensive the- 
ory was elaborated :20., 2l]. Light trapping at the gain in- 
side an RPhC with subsequent release due to the collision 
with another GS have also been predicted and correspond- 
ing memory device was also suggested [32]. In this paper, 
we look at the two separate effects associated with a local- 
ized population inversion inside RPhCs: higher-order soli- 
ton dynamics in the presence of such gain defects and gen- 
eration of the gain defect itself. This paper is subdivided as 
follows. In Sect. [2] the bi-directional Maxwell-Bloch model 
is formulated and analytical results are given for the exact 
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GS interacting with the localized gain inside the RPhC. 
Sect. |3] contains numerical results for the interaction of 
a GS train with this inversion length. Sect. [4] discusses 
possible experimental scenarios and presents conclusions. 



2 Analytical Formalism 

In order to explore the basic features of interaction of res- 
onant gap solitons in the presence of a localized gain or 
defect inside the photonic crystal we consider the simplest 
cases where the defect is created in the middle of the RPC 
using an incoherent pump with square hump-like spatial 
profile 



2.1 Coupled-Mode Equations 

The ID resonant photonic crystal is assumed to comprise 
of periodically positioned thin layers of two-level atoms. 
The electric field that appears in the Maxwell's equations 
and couples the two levels is taken in the form of a super- 
position of the forward, E + , and backward-propagating 
wave, E~ , as follows: 

E(x,t) = ^[E+(x,t)e l( - kx -^ +E-{x,t)e- l[kx+ult) ]{l) 

The direction of propagation is normal to the layers. The 
incident radiation has a wavelength A which exactly sat- 
isfies the Bragg condition d — A, where d is the RPhC 
period. We also assume that the field described by Eq. ([TJ 
is tuned into the exact resonance with atomic absorption 
of the two-level system (which can represent atoms, ions, 
or excitons). In this case the dynamics of a two- level sys- 
tem, can be described by the generalized Bloch equations 
as follows jSlirapmi^E^PS]: 

P t {x h t) = n(x h t) [n+{ Xl ,t)e lkx > + n-( Xi , t)e~ ikXi ] ,(2) 
n(xi,t) = -Rc{p*(xi,t)[f2 + (x t ,t)e ikx ' 

+ Q-{x h t)e- ik *<]), (3) 

where [2 ± — 2t c \jlE ± /fi with E ± from Eq. |l]) and t c = 
(ZireTj/ZcpX 2 ) 1 / 2 i s the cooperative time that is the mean 
photon lifetime in the structure, Tj is the excited state life- 
time, e is the dielectric constant of the medium, p is the 
density of the resonant atoms, p, is the transition matrix 
element, Xi is the position the i-th resonant layer, c is the 
speed of light in vacuum, and the subscript t indicates the 
corresponding derivative. In turn, the resonant polariza- 
tion P defines the dynamics of the optical field D ± inside 
the medium via the set of the two-mode equations: 

f2±(x,t)±f2±(x,t) =J2^ ikXi P(xi,t)6{x-Xi), (4) 

i 

which after averaging over the structure period forms the 
closed set with Eq. ^ and Eq. pj. This is shown here as: 



nf ± a 



± 



P, 



P t = n(Q + + Q~) 
n t = -Re{P*(Q + 



(5) 
(6) 
(7) 



where the subscripts x and t stand for the corresponding 
partial derivative along the dimensionless time and prop- 
agation coordinate t — t/r c and x — x/(ct c ). As input 
values, we use the following boundary and initial condi- 
tions: 



Qf{x = Q,t) = Y^^i sech[0 - U)/t p ], 



S7~(x = L,t) 
Qf(x,t = Q) 
P(x,t = Q) 

n(x, t = 0) 



1, \x-L/2\<l/2, 
-1, \x-L/2\>l/2, 



(8) 

(9) 
(10) 
(11) 
(12) 



where Tp is the normalized duration of the incident pulse. 
Eqs. ((8|-([l2]) together with Eqs. fully define the 

system under considertion. 



2.2 Lagrangian and mechanical analogy 

The invariable quantity of Eqs. (|5|-((7| is accessible via a 
transformation to the Bloch angle ff[x,t), the sum U = 
f2 + + Q~ and the difference V = fi + — Q~ fields, which 
leads to the following set of equations: 



Vf 



- V x = -2 sin 6, 
U x = 0, 
6 t = U. 



(13) 
(14) 
(15) 



An invariable quantity r(x) associated with t his set of 
equations can be obtained by substituting Eq. ( 15 ) into 
Eq. (14 1 and and results in 



r(x) = v(x,t) + e x (x,t). 



(16) 



This conservation parameter is defined by the initial 
conditions and does not change as the pulse propagates 
through the RPhC. However, it does describe the behavior 
of the GS in the sense that setting it equal to zero provides 
a GS whereas, a non-zero value of r(x) leads to a lost of 
stability of the GS and can result in trapping of the GS 
inside the RPhC. It is readily seen that the substitution of 
the invariant function described by Eq. (16 1 into Eq. (13 1 



transforms it into the perturbed Sine-Gordon equation for 
the Bloch angle 8 shown as: 



-9,, =2 sin* 



(17) 



In order to elucidate the interaction between the exact 
GS and defect mode that is due to the localized gain, we 
apply the approach developed earlier in |23j where the 
coordinate of the GS center-of-mass ((t) can be derived 
from a corresponding Lagrangian density: 



L= l -e 2 t - l -(e x -r) 2 -(i-cos0), 

with corresponding density of the Hamiltonian 



(18) 



H 



1 



1 



1 



+ ^i-re x + -r 2 + (i-cos#), (19) 
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where the first four terms on the right-hand side define the 
energy density (U 2 + V 2 )/2 of the forward and backward 
waves inside the RPhC. 

Assume now that the solution of Eq. (17) can be ap- 



proximated by that of the non-perturbed sine- Gordon equa- 
tion: 



, t) = 4 arctan exp 



V~l - u 2 (t) 



(20) 



where and u(t) are the time-dependent coordinate 
of the GS center and velocity, correspondingly. Since the 
total energy of the localized solution is an invariable quan- 
tity, the energy of the interaction of the soliton of Eq. ( 20 1 



with perturbation created by the localized gain, is defined 
by the overlap integral [20,21,22125]: 



<2> = 



1 



+00 



dxF x t 



(21) 



which in turn describes the GS propagation across the 
RPhC with a length of the gain inside as a classical motion 
of the particle with unit mass, 



Cm = 



(22) 



Our initial conditions of the gain localized in the mid- 
dle of the RPhC, namely: 



9(x,t = 0) 



< x < (L-£)/2, 
{L-l)/2 <x< (L + £)/2, (23) 
< (L + l)j2 < x < L, 



along with the absence of the propagating waves in the 
beginning, generates the following potential: 



(p 



-sech [x - (L + i) /2] + sech [x-{L- I) /2] . (24) 



The profile of this potential and corresponding phase tra- 
jectories are depicted in Fig. [I] it is readily seen that if 
the velocity of the particle (pulse) is considerably above 
the potential zero level, it experiences some minor varia- 
tions while traversing the potential, or gain hump in our 
case. Lowering the velocity shifts the phase trajectory into 
vicinity of the saddle point (not shown in Fig. lb) and the 
potential becomes repelling that implies a blockage of the 
GS free transport across the RPhC. Next, if the velocity 
is small enough, the particle experiences aperiodic oscil- 
lations inside the potential well on the left. That is, the 
slower GS is trapped by the gain and oscillates around the 
gain center, it is unsteady but stable. It is worth noticing 
that the energy gain of the GS that is due to the gain of 
Eq. ( 23 ) along with a modification of this gain by the GS 



can change the dynamics of the interaction to a very large 
extend. 



3 Calculated Behavior and Discussion 

The analysis of the previous Section is now used to ex- 
plore the behavior of gap solitons in the presence of a 
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Fig. 1. Profile of inversion n (dotted line) and corresponding 
potential $ (solid line) created in the RPhC for the GS of 
Eq. ( 20 1 (a) along with corresponding phase portrait (b); notice 
existence of turning point, separatrix trajectory, and closed 
aperiodic cycles. 
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Fig. 2. Localization and storage of a single 2n pulse inside 
RPhC on an inversion induced defect; Q\ = l,ti = 50, O2 = 
0,1/ L — 0.02, where I is the length of the inversion and L 
is the RPhC length; time and distance are normalized to the 
coherence time and length, correspondingly. 
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Fig. 3. Release of a stored optical pulse by out-of-phase 
readout 2tv pulse and their merger after (a) 1.5 oscillations 
around the defect (delay time ti — ti = 350) and (b) 2.5 os- 
cillations around the defect defect (delay time t 2 — t / = 400), 
i?i = 1 = — J?2, II L = 0.02; (c) breathing structure of the 
merger pulse — inversion scan at x = 1100 for (b). 



Fig. 4. Merger of the readout pulse with the stored one and 
their pinning by the defect in the form of quasi-symmetric 
bound state; time t% — tj = 475. In this case pulses collide 
on the defect after 5 oscillations of the leading pulse around 
the defect. 



defect mode that is due to the total population inversion 
at the centre of the sample. The class of interactions we 
are concerned with here are those that produce a light 
stopping and release operation which can be exploited in 
order to achieve a memory cell with high-speed switching 
at a moderate consumption of the laser power. 

We have previously shown that it is possible to trap 
a gap soliton on the potential of Fig. 1 [19jl21j. Figure 2 
shows a contour plot in the (a;, i)-plane, depicting the local 
inversion and the propagation and subsequent trapping of 
a 2ir pulse. We call the GS that is stopped and oscillat- 
ing at such defect an information bit. It is important to 
note that only such stopped solitons can be used to switch 
subsequent pulses providing either true memory readout, 
higher-density storage, or limiting action so providing a 
secure readout. This trapping represents a basic optical 
memory whereby the GS is stored in the RPhC. However, 
in order to demonstrate its usefulness we need to inves- 
tigate the read-out or release of information from such a 
memory device. 

In Figs. 3a and 3b we show that a second 2ir gap soli- 
ton can be used to release the stored pulse and effectively 
read out the optical memory. The soliton readout hap- 
pens owing to the collision of the out-of-phase pulse with 
the trapping area. Our simulations show that the soliton 
can be released by the readout pulse as long as the read- 
out GS is timed to arrive when the trapped soliton is at 
its maximum deviation towards the front side of the de- 
fect. As the readout soliton arrives it interacts with the 
trapped pulse and is dragged into the defect. The trapped 
and readout solitons cross the center of the defect together 
and have sufficient energy to overcome the potential and 
are released. Naturally, the released pulse is a higher or- 
der breather (see Fig. 3c) and it propagates with a higher 
velocity through the resonant photonic crystal. 
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From this argument it is clear that if the readout GS 
is not synchronized with the trapped soliton then they 
can combine at the defect and can both be trapped. We 
illustrate this effect in Fig. 4, where the out-of-phase soli- 
ton used to trigger the read-out arrives while the trapped 
one is at its turning point on the opposite side of the de- 
fect. The resulting interaction causes the solitons to arrive 
at the center of the inversion region from opposite sides. 
The result is that the solitons combine to form a trapped 
mode at the inversion. This corresponds to a higher order 
trapped mode as can be seen from the higher frequency 
structure on the soliton envelope. In order to implement 
a practical optical memory it is essential that the solitons 
be timed correctly. This implies that the strength of the 
inversion in the photonic crystal should be such that the 
period of oscillation of the trapped mode is the same as 
the bit rate of the data stream. If this condition is not 
satisfied, the merging breather-like state turns out to be 
unstable, dissociates into a soliton trapped to the inver- 
sion and the read-out pulse rejected backward - a read-out 
failure. This observation (which is not shown in Figures) is 
extremely important from the information security stand- 
point. That is, since the oscillation period and amplitude 
determine the open time of the memory cell a pulse with 
the wrong phase cannot take away the information bit. 

Since the process of the memory read-out is due to 
the energy imbalance between interacting solitons and the 
trapping defect, it is natural to assume that trapping and 
release of a gap soliton may also be controlled by means of 
its interaction with an in-phase, and more intense pulse. 
As a relevant example, Fig. 5 shows the results of numeri- 
cal integration of the Maxwell-Bloch equations, Eqs. ([5])- 
([7| for the case of interaction of two solitons with different 
input amplitudes. Above a certain threshold, the trailing 
GS with larger energy permits the escape of the GS from 
the trapping area (Fig. 5a). This is due to the repulsive 
action of the read-out pulse that is stronger than the trap- 
ping action provided by the defect. What might be inter- 
esting from a potential device feasibility stand-point is 
that the released soliton emerges the RPhC in the back- 
ward direction. Again, this happens only for a certain bit 
rate of the data stream. The readout failure is demon- 
strated in Fig. 5b and 5c. Here, the speed of the released 
light is so low that is practically creates another defect 
inside the crystal so rendering the original one available 
for the information storage. 

The other important finding that follows from Fig. 5 is 
the existence of a stabilization of a higher-order solution to 
the perturbed sine-Gordon equation, Eq. (17), due to the 



presence of a local perturbation. This is in sharp contrast 
to the well-known fact of the breakup of the 27rn-pulse 
(with n being an integer) of the self-induced transparency 
into a train of isolated 2ir pulses with different amplitudes, 
speeds, and durations |2"I] . 

We conclude this Section with indication of an alterna- 
tive way of stopping light inside a RPhC. Contrary to the 
initial conditions of Eqs. (|8j-(12| let us now assume that 
the gain (or total inversion of the atoms) is not squeezed 
in the middle of the sample but is distributed from the 
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Fig. 5. Release of a stored optical pulse by in-phase intense 
pulse (a) and readout failure (b); fli — 2/?2/3,t2 — ti = 250 
(a), ti-ti = 275 (b). 



input facet at certain distance inside. Using initial con- 
ditions pertinent to spontaneous decay [2"S"lErJ] . also leads 
to interesting results. This is summarized in Fig. 6. The 
localized field oscillating around the physical boundary 
between the pumped and unpumped parts of the RPhC 
is clearly visible. 



4 Conclusions 

The results presented in this paper are based on studying 
gap solitons in resonance photonic crystals. However, the 
sine-Gordon equation can also be used to describe a broad 
and diverse range of physical systems; for example, the dis- 
sociation of the breather solution of the sine-Gordon equa- 
tion that is caused by a perturbation of a long Josephson 
junction ^ I2 ^ [29 | [3g i [3l j l32] . The problem of the stability 
of the perturbed sine-Gordon equation also results in the 
ratchet effect [331I34] . the high-frequency field interaction 
with electron plasmas in a semiconductor super lattice |35j , 
and the spatial localization in a nonlinear chain 36J. We 
would therefore expect that similar effects of energy stor- 
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Fig. 6. Spontaneous decay and oscillating defect formation 
at the interface between pumped and ground-state parts of the 
RPhC; the normalized length of the pumped section is 500 (a), 
600 (b), and 800 (c). 



age and release by a defect will also be present in these 
systems. 

In conclusion, we have demonstrated that the oscillat- 
ing immobile gap solitons created by the presence of a 
localized region of the population inversion (or gain) in- 
side a resonance photonic crystal can be manipulated by a 
proper choice of bit rate, phase and amplitude modulation. 
Developing this idea, we are able to obtain qualitatively 
different regimes of the RPhC operation. A noteworthy 
observation is that both the delay time and amplitude 
difference must exceed a certain level to ensure effective 
control over the entire soliton dynamics and its speed, in 
particular. A modification of the defect that accomplishes 
a soliton slowing-down and trapping can make the dynam- 
ics of soliton trains in the resonance photonic crystal with 
defects even more interesting. This comprises a subject of 
our future work. 
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